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The vibrational dynamics of protein folding is analyzed in the framework of Tsallis 
thermostatistics. The generalized partition functions, internal energies, free energies 
and temperature factor (or Debye-Waller factor) are calculated. It has also been 
observed that the temperature factor is dependent on the non-extensive parameter 
q which behaves like a scale parameter in the harmonic oscillator model. As q — > 1, 
we also show that these approximations agree with the result of Gaussian network 
model. 

Keywords: Tsallis thermostatistics, Harmonic Oscillator, Protein Folding, Tem- 
perature Factor 

PACS numbers: 05.20.-y,87.10.te,05.70.-a 



INTRODUCTION 



The statistical mechanics paradigm based on Boltzmann-Gibbs entropy which is great 
success looks after be incapable to deal with many interesting physical systems[l|]. In 
1988, Tsallis advanced a nonextensive generalization of Boltzamann-Gibbs entropic mea- 
sure. There are many applications for Tsallis formalism such as the specific heat of the 
harmonic oscillator 2], one- dimensional Ising model jjj, the Boltzmann H-theorem {4], the 



Ehrenfest theorem 



5j, quantum statistics (6j, paramagnetic systems Cercle maps 
Henon map [9], Haldane exclusion statistics lOj, q-expectation value 11], quantum me 



chanical treatment of constraints 



121 ]. In this paper, we will focus on thermodynamics and 
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vibration properties of proteins in the framework of a generalized thermostatistics. Before 
proceeding with generalized gaussian network model, we described fundamental definitions 
of Tsallis thermostatistics in Section II. We then study generalized gaussian network model 
and result in section III. Results will be discussed in section IV. 

II. FUNDAMENTAL DEFINITIONS OF TSALLIS THERMOSTATISTICS 

The generalization which has been successfully nonextensive statistical mechanics is based 



on the following expression 13j 



1 _ Tr(o q ) 

s q = k l q _[ p) . (1) 

In the Tsallis statistics, the nonextensive canonical distribution is given by 

Pq = ^[l-(l-q)(3H} 1/{1 - g \ (2) 

where 

N 

Z q = / J] d Pndx n [l - (1 - q)(3H} 1/(1 ~ q \ (3) 

J n=l 

is the partition function with q G R and f3 = 1/T (with k = 1). Eq.(2) is obtained 



13, 



14j | . The internal energy 



maximizing the Tsallis entropy using lagrange multiply methods 
is described by 

N 

U q = Hdp n dx n p«H (4) 

J n=l 

where 

N 

/ J]dp n da; fl p« = l. (5) 

n=l 

In Free energy notation, Eq.(4) appears as 

U q = F q - TdFjdT, (6) 
where F q is free energy is described as 

F q = {Z*- 1 -1)/[(1 -q)0\. (7) 
The expectation of an observable A is 

N 

A q = (A) q = / J] d Pndx n p q - (8) 

J n=l 
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III. GENERALIZED GAUSSIAN NETWORK MODEL AND RESULTS 



We focus our discussion on the classical Tsallis statistics by considering protein folding 
and harmonic approximation. Consider the following Hamiltonian for Gaussian Network 
Model(GNM) based on harmonic approximation 15|] 



H 



N „2 

E 



i=l 



p- + ^-ARfLARi 
2m 2 



(9) 



where the first term is the kinetic energy of the system and 7 is the strength of the springs. 
In this model, springs are assumed as homogeneous. Ri and ARi is also defined as the 
equilibrium position and the displacement with respect to Ri of the i-th C a atoms. The 
model is eventually described by the contact matrix L with entries: Ly = 1 if the distance 
|Ri — Rj| between two C a 's, in the native conformation, is below the cutoff Rq, while is 



otherwise. From Eq.(3), partition function for the Gaussian Network model 15] which has 
above hamiltonian is given by 

N 



Z q = f f[d K dA^[l-(l- g )/5 

' N 

E 



jv 2 

+ -ARfLARi 
2m 2 

i=l 



(10) 



we introduce the variables L = V T XV, x = VAR and thus AR T LAR = x T \x if A is 
eigenvalues of the L and V is eigenvectors of L. We can rewrite partition function as: 

N 



Z q = I JJdptdaritl-tl-g)/? 



N 



1 1/(1-9) 



X 



,i=l 



(11) 



To calculate this integral, we definition new variables: = [(1 — q)j\iP/2] 1 ^ 2 Xi and Un+i — 
[(1 — q)f3/(2m)] 1 / 2 pi [16l |. where i = 1,2, 3, N. Substituting these variables into Z g gives 



N r 



n 



,m 



x 



2N / 2N 

n*» i-E 



n=l 



k=l 




V(l-9) 



(12) 
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By using hyperspherical coordinates with u = (X^=i Un) 1 ^ 2 an d calculating the integral over 
the angular variables [la], we obtain 

N 



n 

.71=1 



1/2 



(l-g)A n /r 7 



n 



2N 



x / duu N ~ l (l-u) 1/{1 - q) 



(13) 



By using definition of the solid angle 
of Euler beta function [3], we obtained that 



171 ] Q2N = 2n N /T(N) and using integral representation 



N 

n 

n=l 

2-g 
1-g 



27T / m \l 



(l-g)A n /T 7 

1 JV 



/2 



rf^ + i + iv 



AT. 



n 



n=l 



2tt 



,m ly/2 



|_(l-g)A n /3 v 7 



(14) 



where (a) n = a(a + l)(a + 2). ..(a + 71—1) 
energy of the system is given by 



^[a-iy l 1S ^e Pochhammer symbol [18j. Free 



1 



F q - --\n q Zq, 



(15) 



and 



2-g 

l-g 



AT. 



-1 V 



n 

n=l 



2tt 



'™Nl/2 



|_(l-g)A n /? v 7 



(16) 



and then 



2-g 

l-g 



A?J 



27T m j 



/2 



L(l-g)/5 v 7 



A' 



(det {L )) 



UN 1/2 



(17) 



From Eq.(6), the generalized internal energy can be calculated as 



' /3 ' /3 



27T ,m. 



1/2 



;i-g)/5 v 7 



(det(L- 1 )) 



1NN1/2 



2-g 
i-g 



A. 



-1 



1-9 



(18) 



Hence, we obtained thermodynamic properties of the generalized Gaussian network model. 
On the other hand, the most properties of Gaussian network model is called temperature 
factor. The crystallographic temperature factors or B-factors, which is defined as intrinsic 
fluctuations of the atoms in crystal, is directly related to this fluctuations. The X-ray 
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crystallographic temperature factors are used by the mean square fluctuation of C a atoms 



around their native positions 



151 ] . The generalized temperature factor is expressed as 



B iq {T) 



57T 



with (•) indicating the thermal average. Where {xiXj) q is defined by 

(xiXj) q = -^g J DxDpXiXj[l - (1 - q)pH}^, 



or 



(XiXj^q q / JDxDpX'iXj 



N 



1=1 



1 

1-9 



(19) 



(20) 



(21) 



In the GGNM (generalized gaussian network model) approximation, this average is easily 
carried out, because amounts to a nonGaussian integration, and B-factors can be expressed 
in terms of the diagonal part of the inverse of the matrix 



(Xi ■ Xj)c 



3k B T, 



7 L l + (l-g)(iV+l)- 



(22) 



Substituting Eq. 20 into Eq.17 gives 



7 L l + (l-g)(iV+l) 
If q — > 1, it is calculated that same result of the gaussian network model as[l9j 

8n 2 k B T 



(23) 



Bi(T) 
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(24) 



For the numerical results, we choice a reasonable cutoff distance including all residues pairs 
within a first interaction shell is 7.0 A. We display in Figfj] temperature factors for T4 
lysozyme (or 3LZM pdb) as a function of residue index at q=0.99. Bold and dashed curves 
depict the experimental and theoretical results, respectively. The theoretical curve is nor- 
malized by taking 7(A~ 2 ) = 2.06&bT, therefore as to match the area enclosed by the two 
curves. The agreement between the theoretical and experimental curves is remarkable. In 
Fig 121 we also display to behavior cross correlation between i. residues and j. residues of all 
modes at q=0.99. 
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FIG. 1: Temperature factors for T4 lysozyme 

IV. CONCLUSION 

In this study, we obtained temperature factor or Beta factor using by Tsallis statistical 
mechanics. In this framework, by the approximate scheme which we have used here, closed 
and analytical expressions for the generalized gaussian network model have been derived 
for harmonic oscillator. In addition to this, the temperature factor is q dependent. Since 
the temperature factor must be positive, q must be less than unity. This also shows that q 
must be less than unity for harmonic oscillator system. As can be seen in FigJT], at q=0.99, 
the beta factors which are obtained from experimental and numerical calculation are in 
agreement with each other. The more the q parameter decreases, the more the temperature 
factor disagrees with experimental data. Therefore, for q values in the lower vicinity of 
1, the experimental data are in agreement with numerical calculations. In this system, as 
q — > 1, it is observed that, both of the above expressions reduce to the results obtained in 
the standard Boltzmann-Gibbs statistics. 
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FIG. 2: Cross correlation between pair residues for T4 lysozyme 
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